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ABSTRACT 

o : 
o . 

CN ' We report on optical spectroscopy of the eclipsing Halo LMXB 2S 0921- 

630/V395 Car, that reveals the absorption line radial velocity curve of the KOlll 

XJ^ • secondary star with a semi-amplitude i^'2=92.89± 3.84 kms~^ , a systemic ve- 

O ■ locity 7=34.9±3.3 kms'^and an orbital period Porb of 9.0035±0.0029 day (l-a). 

^ ■ Given the quality of the data we find no evidence for the effects of X-ray ir- 

^ ■ radiation. Using the previously determined rotational broadening of the mass 

^ ' donor, and applying conservative limits on the orbital inclination, we constrain 

^ ■ the compact object mass to be 2.0-4.3 MQ^l-a), ruling out a canonical neutron 

Q\ • star at the 99% level. Since the nature of the compact object is unclear, this 

^ ' mass range implies that the compact object is either a low-mass black hole with 

O ■ a mass slightly higher than the maximum neutron star mass (2.9 Mq) or a mas- 

• sive neutron star. If the compact object in 2S 0921-630/V395 Car is a black hole, 

it confirms the prediction of the existence of low-mass black holes, while if the 

2 ■ object is a massive neutron star its high mass severely constrains the equation of 

^ ■ state of nuclear matter, 
cd ■ 

> ' 

■ Subject headings: accretion: accretion discs - binaries: close stars: individual 

^' (2S0921-630/V395Car) 
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1. Introduction 

A knowledge of the neutron star mass distribution provides a fundamental test of the- 
ories of the equation of state of nuclear matter, the applicability of General Relativity as 
the correct theory of gravity and significant information on the evolutionary history of the 
progenitor stars. To date, only studies of millisecond radio pulsars have provided accurate 
mass determinations of neutron stars, reflecting their mass at formation; 1.35±0.04 Mq (van 
Kerkwijk 2001). Although dynamical neutron star masses can also be obtained from accret- 
ing X-ray pulsars, there are large uncertainties. In high-mass X-ray binaries (HMXBs) there 
are uncertainties due to non-Keplerian perturbations in the radial velocity curves caused by 
effects such as stellar wind contamination and X-ray heating. In low-mass X-ray binaries 
(LMXBs), the situation is worse because the intense X-ray irradiation usually suppresses the 
light from the donor (Charles & Coe 2003). It is only in a few exceptional cases where the 
companion is evolved (and hence more luminous) or during X-ray off-states, when dynamical 
information can be extracted about the nature of the compact object. 

The SAS-3 X-ray source, 2S 0921-630, was identified with a ~16™ blue star (Li et al. 
1978), V395 Car, whose optical spectrum is dominated by HeII4686A and Balmer emission 
(Branduardi-Raymont et al. 1983), characteristics of LMXBs, with occasional superimposed 
late- type stellar absorption features visible (Thorstensen & Charles 1979). EXOSAT obser- 
vations showed a broad (>1 day), shallow X-ray eclipse, during which the spectrum softened 
(Mason et al. 1987), allowing an estimate of the size of the accretion-disc corona (ADC). 
This was comparable to the eclipsing companion star, thus accounting for the partial X-ray 
eclipse, and with an inclination of i=70°-90°. Only a handful of ADC sources are known in 
which the compact object is permanently obscured from our line-of-sight by the accretion 
disc, and requires the observed X-rays to be scattered in a hot corona (White et al. 1995). 
By implication, the intrinsic X-ray luminosity is much higher. 

Limited optical photometry and emission line spectroscopy (Cowley et al. 1982; Branduardi- 
Raymont et al. 1983) shows that 2S 0921-630/V395 Car has a long orbital period of 9.02 day 
and deep (~1 mag) dips where the {B — V) colour reddens by up to 0.4 mags, which has 
been interpreted as the echpse of the disc by the late-type secondary (Chevaher & Ilovaisky 
1982). High-resolution optical spectroscopy revealed the donor star to be of spectral type 
KOlll with a rotational velocity of vsin i=65±9kms~^, contributing ~25% to the observed 
flux at 6500A (Shahbaz et al. 1999). 2S 092 1-630/ V395 Car is one of those rare LMXBs in 
which the secondary is visible despite the presence of a luminous disc. 

There has been no detection of type I X-ray bursts, so the nature of the compact 
object is unclear. In this letter we present the results of an intensive campaign to measure 
the dynamical mass of the compact object in 2S 092 1-630/ V395 Car. We determine the 
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secondary star's radial velocity curve, investigate the possible effects of X-ray irradiation 
and obtain firm constraints on the compact object mass. 

2. Observations and Data Reduction 

Time resolved spectroscopic observations of 2S 0921-630/V395 Car were obtained on the 
1.9-m telescope at the South African Astronomical Observatory (SAAO) during 2003 March 
30 to April 18. Grating # 5 (1200 lines/mm) was used centered at 4900A and covering a 
wavelength range of 4529-5377Awith a dispersion of 0.49 A/pixel~^. A slit width of 2.0"gave 
a spectral resolution of 1.1 A (=64 kms~^ at 5200A) as measured from the Cu-Arc arc lines. 
In general the conditions were good, with the seeing varying between 1 "-3 A total of 
83 spectra of 2S 0921-630/V395 Car were obtained with an exposure time of 1800 s. Each 
spectrum was bracketed by an observation of the internal Cu-Ar arc lamp at the position of 
the star. Template field stars of a variety of spectral types were also observed. 

The images were fiat-ficldcd using observations of a tungsten calibration lamp, and 
a constant bias level was subtracted from all the images. Extraction of the 1-D spectra 
from the images were performed using optimal extraction (Horne 1986) and calibration of 
the wavelength scale was achieved using 4th order polynomial fits to the position measured 
arc-line positions, which gave an rms scatter of ~0.02 A. 

In addition, 34 spectra were also obtained at the Very Large Telescope (VLT; Paranal, 
Chile), New Technology Telescope (NTT; La Silla, Chile) the Anglo Austrahan Telescope 
(AAT; Siding Springs, Austraha) and at the Magellan telescope (Las Campanas, Chile) as 
part of a long-term project on X-ray binaries. We used the F0RS2 Spectrograph attached to 
the 8.2m Yepun Telescope (VLT) on the nights of 2003 May 18 and June 22-23. A total of 9 
spectra with exposure times in the range 300 to 900 s were obtained. The R1400V holographic 
grating was used in combination with a 0.7 "slit resulting in a wavelength coverage of 4514- 
5815A and a resolution of 70 km s^^. Instrumental flexure was monitored by cross-correlating 
the sky spectra and was found to be very small, always within 6 km s~^ . The Magellan spectra 
were obtained on 2003 Dec 12-16 using the Baade telescope at Las Campanas using the 
IMACS imaging spectrograph in long-camera mode. A total of 11 spectra with an exposure 
time of 300s were obtained. The 600 grating together with a 0.7 "long slit provided spectra 
covermg 3720-6830A with a dispersion of 0.75A/pixer^ The AAT spectra were obtained 
using the RGO spectrograph on the nights of 2002 June 6 to 11. A total of 10 spectra with an 
exposure time of 1800 s were obtained. The R1200B grating centered at 4350A together with 
a 1.0 "slit provided spectra covering 3500-5250A with a resolution of 70kms^^. The total of 
4 NTT spectra were obtained with EMMI and grating # 6 on the nights of 2002 June 8 to 
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10 exposure times of 600 to 900 s. The spectra cover the spectral range 4400-5150A with a 
resolution of 75kms~^. The AAT and NTT data reduction details are given in Casares et 
al. (2003). 

3. The radial velocity curve and spectrum 

To increase the signal-to-noise of the spectra prior to cross-correlation, the individual 
spectra were variance averaged into nighly means or groups of ~3 hrs, depending on the qual- 
ity of the individual spectra. A total of 31 absorption line radial velocities were measured by 
cross-correlation (Tonry & Davis 1979) with a template star. Prior to cross-correlation, the 
spectra were interpolated onto a constant velocity scale (32 kms~^ pixel"^) and normalized 
by fitting a spline function to the continuum. The region 5100-5300A common to all the spec- 
tra, primarily containing the Mgl (5167.3A, 5172. 7A and 5183. 6A) absorption blend was used 
in the analysis. The template star's radial velocity, determined using the position of the Fel 
4957. 597A absorption hue, was then added to the radial velocities of 2S 092 1-630/ V395 Car. 
The resulting radial velocity data (see Figure 1) were then fitted with a circular orbit (si- 
nusoidal function) using the method of least-squares. Using spectral type template stars 
in the range G6-K4III we obtain K2 values in the range 92-99 kms~^(see Table 1). The 
secondary star's spectral type has previously been determined from high-resolution spec- 
troscopy to be KOlII (Shahbaz et al. 1999), and our work confirms this. Fitting the radial 
velocities with a KOlll template star, we find a minimum reduced of 4.3 with the best 
fit parameters 7 = 34.9±3.3 kms'^ K2 = 92.89±3.84 kms"^ Porb=9-0035±0.0029 day and 
To=HJD2453099.51±0.08, where Porb is the orbital period, Tq is time at phase 0.0 defined 
as inferior conjunction of the secondary star, 7 is the systemic velocity and K2 is the ra- 
dial velocity semi-amplitude (l-cr errors are quoted with the error bars rescaled so that the 
reduced of the fit is 1.) Our ephemeris is consistent with that determined by Mason et 
al. 1987. We also computed a periodogram of the data to investigate the significance of 
other periods. No significant peaks (at >99.99 % level) were present other than the one at 
9.0035 day. Finally, to check if the 9.0035 day period is affected by an alias, we computed the 
window function; no significant peaks were present at 9.0035 day. Fitting the radial velocity 
data with an eccentric orbit did not yield a better fit; the eccentric fit is only significant at 
the 38% level. 

Using a distance of ~10kpc Cowley et al. (1982) 2S 0921-630/V395 Car lies 1.6 kpc 
below the galactic plane and so belongs to the Halo population. Note the systemic velocity 
is not consistent with the radial velocity due to Galactic differential rotation, ~60 kms"^ 
at a distance of ~10kpc (Dehnen & Binney 1998). However, what is interesting is that the 
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magnitude of the systemic velocity of 2S 0921-630/V395 Car seems to be about a factor of 
~4 lower compared to other neutron star Halo systems (Casares et al. 1998, 2002; Torres et 
al. 2002). 

In Figure 2 we show the variance- weighted Doppler averaged spectrum of V395 Car. 
The strongest features are the H, He emission lines, which appear single peaked, even at 
the resolution of these data. Weak absorption features of Mg I 5175A and the Ca/Fe blend 
at 5269A allow us to determine the fraction of light from the secondary star, by optimally 
subtracting a scaled KOlII template spectrum from V395 Car (Marsh et al. 1994). As the 
rotational broadening of the donor (65±9 kms~^; Shahbaz ct al. 1999) is similar to our 
spectral resolution (64 kms^^), we do not broaden the template before subtraction. We find 
that the secondary star contributes ~ 10% to the observed flux at 5200A (see Table 1). 

4. Irradiation of the secondairy star - the mass of the compact object 

The observed X-rays are reflected into our line of sight by scattering. Although the 
observed X-ray luminosity is low Lx = 2 x 10^^ erg s^^ (Kallman ct al. 2003), it is still 
interesting to explore the possible effects that X-rays may have on the secondary star's 
atmosphere and hence radial velocity. We investigate these effects using the X-ray binary 
model described by PhiUips et al. (1999); Shahbaz et al. (2000, 2003). We flrst phase-folded 
the radial velocities using the orbital ephemeris derived in section 3 and binned into 31 orbital 
phase bins. We then flt the resulting radial velocity curve with the X-ray binary model. We 
considered two extreme cases for the X-ray irradiation: case (a) where the X-rays produced 
near the compact object irradiate the whole inner face of the secondary star and case (b) 
where the X-ray's do not affect the secondary star. Note that invoking a flared accretion 
disc will have intermediate solutions. The variable model parameters are the compact object 
mass Ml, secondary star mass M2, X-ray luminosity Lx and the inclination i. To determine 
which elements on the star contribute to the absorption line radial velocity we use the factor 
fx, which is the fraction of the external radiation flux that exceeds the unperturbed flux. 
For a given i and combination of Mi and M2 we flt the data and use the constraints imposed 
by the measured value for vsin i (64±9kms~^; Shahbaz et al. 1999) to produce plots in the 
M1-M2 plane (Figure 3) 

Given the evidence for optical and X-ray partial eclipse, i must he in the range 70°-90°. 

Therefore we consider cases (a) and (b) for i=70° and more importantly i=90°, because it 
gives a firm lower limit to the compact object's mass. For case (a) we fix Lx = 2 x 10^^ ergs~^ 
and find that the best fit (x^ = 122) is obtained when all the elements on the inner face 
contribute to the radial velocity, i.e. fx > 60%, i.e. X-ray heating is not preferred. Note, 
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however, that the radial velocity curve will not be sinusoidal because there will still be 
some irradiated elements that contribute to the radial velocity. Even if we increase to 
Lx = 10'^*^ erg s"^ as could be the case for an LMXB, the fits are significantly worse (x^ = 
190) and the best fit still occurs when all the inner face contributes to the radial velocity. 
For case (b) we set Lx = 0, implying no X-ray heating and find that it gives a fit with 
similar quality (x^ = 119). We conclude that the best fit to the radial velocity curve is the 
case where the X-rays do not have an observable effect given the quality of the data. For 
our range in i our fits for case (a) and (b) give 2.0 < Mi < 4.3 Mq (l-cr), which means 
we can rule out a canonical 1.4 M© neutron star at the 99% level. Our solutions for the 
compact object mass suggests either a massive neutron star or a low-mass black hole, since 
the maximum mass allowed for a neutron star is 2.9 M© (Kalogera & Baym 1996). 

Since the donor star fills its Roche lobe and is synchronized with the binary motion, we 
use the secondary star's v sin i and K2 to determine the binary mass ratio g'=(M2/Mi)=0.89±0.18 
(Horne et al. 1986). Given the measured masses and orbital period, we use Kepler's Third 
Law to determine the semi-major axis 56. 1< a <72.4R0. Eggleton's (1983) expression 
for the effective radius of the Roche-lobe then determines the radius of the secondary 
10. 4< i?2 <13.4R0, and the temperature inferred from the spectral type (Gray 1992) 
and Stefan-Boltzmann's law determines the luminosity 52. 5< L2 <87.1L0(l-(j ranges are 
quoted) . 

5. Discussion 

5.1. The evolutionciry status of the system 

The position of the secondary star of 2S 092 1-630/ V395 Car in an HR diagram corre- 
sponds to a normal star that has crossed the Hertzsprung gap and now lies on the Hayashi 
line. The evolution of the binary is dominated by the evolution of the evolved secondary star. 
Since such a star no longer burns hydrogen in its core, the mass transfer is early massive 
Case B. For a binary with a secondary star near the onset of such mass transfer, one requires 
q <1 throughout its history (King & Ritter 1999). The case for a 1.4 Mq neutron star can be 
ruled out since it requires M2 <1.4 Mq, which is inconsistent with our results; the position 
of a secondary on the HR diagram undergoing early massive Case B mass transfer is close to 
that of a single star of the same mass (Kolb 1998). For a 3 Mq compact object, M2 <3.0 Mq, 
which is consistent with our observations. 

It is interesting to note that the mass ratio we derive lies on the limit of stable/unstabe 
mass transfer g=5/6 (King & Ritter 1999). It is possible that 2S 0921-630/V395 Car is 
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undergoing unstable mass transfer similar to the supersoft sources and other X-ray binaries 
(Kahabka & van den Heuvel 1997). 

5.2. The formation of the compact object 

The distribution of measured neutron stars masses provides fundamental constraints 
on the equation of state of nuclear matter and the observational identification of a black 
hole which is based on the maximum gravitational mass allowed for a neutron star; 2.9 Mq 
(Kalogera & Baym 1996). Studies of radio pulsars show that neutron star masses are clus- 
tered in a narrow range. However, there is some evidence that neutron stars with a mass in 
excess of 1.4 Mq do exist. LMXBs are expected to contain massive (~2 Mq) neutron stars, 
resulting from the accretion of a considerable amount of material over extended (10^ yrs) 
periods of time (Zhang et al. 1985). However, while a neutron star mass of 1.78±0.23 M© 
for the LMXB CygX-2 is reported (Orosz & Kuulkers 1999), Titarchuk & Shaposhnikov 
(2002) have recently proposed a mass of 1.44±0.06 M0 based on the spectral and temporal 
properties of the type I X-ray bursts. The analysis of optical data for the X-ray pulsar 
VelaX-1 suggests that the pulsar has a mass of 1.871^:?^ Mq (Barziv et al. 2002). Proba- 
bly the strongest case for a massive neutron star or low-mass black hole is in the eclipsing 
HMXB 4U 1700-37 which contains a 2.44±0.27Mo compact object (Clark et al. 2002). Such 
masses already test soft nuclear equations of state (Miller et al. 1998a). We obtain a mass of 
2.0-4.3 Mq for the compact object in 2S 0921-630/V395 Car, which lies between the range 
of masses observed for neutron stars and black holes. Models for nuclear equations of state 
which include the effects of three nucleon interactions and realistic models of nuclear forces 
limits the maximum mass of neutron stars to be below 2.5 M0(Akmal et al. 1998). However, 
it should be noted that unconventional forms of matter such as Q-stars do allow the existence 
of extremely massive neutron stars (Miller et al. 1998b) Although our results do not allow 
us to distinguish between a massive neutron star or a low mass black hole, the existence of 
a >2.5 M0neutron star would place strong constraints on high density nuclear matter. 

Theoretical predictions by Fryer & Kalogera (2001) suggest that the death of a 60 Mq 
star in a close binary system can produce anything from a low-mass 1.2 Mq neutron star 
to a 10 M0 black hole, depending on the wind mass-loss rate during the Wolf-Rayet phase. 
As suggested by Brown et al. (1996) the formation of a low- mass black hole with a high 
space velocity is possible via a two stage process involving the formation of a neutron star. 
However, if there is symmetric mass ejection in the supernova, then a kick at the formation 
of a black hole is not required to explain the varied range in the observed space velocities of 
the black hole LMXBs (Nelemans et al. 1999). 
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The unusual mass for the compact object in 2S 0921-630/V395 Car and the fact that 
the nature of the compact object is not clear, conjours up many formation scenarios. If 
it is a neutron star, then its present mass can be explained by either the accumulation of 
matter (1-2 Mq) by a canonical 1.4 Mq neutron star, accreted at a high mass accretion 
rate over a long period of time, or the direct collapse of a massive star forming a massive 
neutron star. If the compact object in 2S 0921-630/V395 Car is a black hole then its low 
mass can be explained by the accretion induced collapse of a massive 2.9 M0 neutron star or 
the formation of a low- mass black hole (Brown et al. 1996). Given the fact that we do not 
observe a peculiar systemic velocity comparable to other neutron-star LMXBs, this suggests 
that we can rule out canonical neutron star formation scenarios where a kick is produced 
during a type II supernova. The most likely scenario for the formation of the compact object 
is the direct formation of a massive neutron star or a low-mass black hole. 
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Table 1: Fits to the radial velocity curve of 2S 092 1-630/ V395 Car (l-cr errors are given). 
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Fig. 1. — The heliocentric radial velocity curve of the secondary star. The crosses show the 
SAAO data and the stars the VLT, NTT, AAT and Magellan data. The solid line shows 
a sinusoidal fit to the data. The data have been folded on the orbital ephemeris and 1.5 
orbital cycles are shown for clarity. 



-13- 



1.4 



0) 

X 



1.2- 



o 



m 
■P 

.iH 

£^1.0 

cS 
■p 



0) 



^^^^ fV- 

Q) 0) no Q) qj "5 

O _ 




KOIII 



0.2 



1 



1 



4600 



4800 5000 
Wavelength (A) 



5200 



Fig. 2. — Top: the variance-weighted Doppler shifted average spectrum of 2S0921- 
630/V395 Car in the rest frame of the secondary star. The most noticeable features are 
the H and He emission hues arising from the accretion disc. Bottom: the KOlll spectral type 
template star HD82565. The spectra have been normahzed and shifted vertically for clarity. 
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Fig. 3. — The effects of X-ray heating on the secondary star's radial velocity curve and hence 
the allowed binary component masses. The left and right plots show the case without and 
with X-ray heating respectively. The solid and dashed contours show the 68% confidence 
level obtained by fitting the radial velocity curve with the X-ray binary model using i=90° 
and 70° respectively (see section 4) The dotted horizontal lines mark 1.4 and 2.9 Mq, 
and represent the canonical and maximum neutron star mass respectively. 



